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Abstract. We give a gauge description of the adiabatic charge pumping in closed 
systems, both in Abelian and non-Abelian processes, and by means of asymptotic 
Wilson loops in a suitable parameter manifold. Our geometric formulation provides 
new insights into this issue, and a very simple algorithm for numerical computations. 
Indeed, as we show first, discretized Berry- Wilczek-Zee holonomies are easy to 
implement. Finally, we study non-Abelian pumping in a solvable four-state model, 
already used in several contexts, to demonstrate the relevance of our approach. 
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1. Introduction 

Geometry and topology play a central role in modern physics [1], obviously in advanced 
areas of field theories, but also in "standard" quantum mechanics. Indeed, since 
Berry's seminal paper [2], an increasing number of quantum phenomena have been 
well understood through the mathematical apparatus of gauge theories [3,4]. Some of 
the most outstanding examples are certainly to be found in quantum transport physics: 
Aharonov-Bohm effect, Thouless pumping, various quantum Hall effects, etc. The 
geometrical and topological features in question arise directly from Hilbert spaces. They 
are no more seen as simple vector spaces, but as bundles based on their own projective 
space [5] or, if the system is parameter-dependent, on a parameter manifold [2,6]. 

The latter construction gives a great insight into physics of adiabatically driven 
quantum systems. Here, we only need to look at a single level and the initial bundle can 
be considerably reduced. For an n-fold degenerate level (n ^ 1), a \]{n) gauge structure 
naturally appears over parameter manifolds. The physical information is then entirely 
carried by Berry [2] (n = 1) or Wilczek-Zee [7] (n > 1) connections. Indeed, they 
induce parallel transports in fibres which are tantamount to the Schrodinger equation 
and determine the adiabatic wavefunction along parameters evolutions (up to an over- 
all dynamical phase factor). Above all, oriented loops in parameter manifolds provide 
parallel transport maps called holonomies. In a standard gauge theoretic language, once 
given a gauge fixing, holonomies are represented by U(n) matrices known as Wilson 
loopsX- Just as representations of physical observables, these matrices transform covari- 
antly with respect to loops' base points. Holonomies are de facto susceptible to enclose 
some "elements of reality" . For example, ongoing research in quantum information tries 
to use them concretely for implementation of quantum logic operations [8,9]. Besides, 
the gauge-covariant curvatures, naturally induced by Berry-Wilczek-Zee connections, 
may for their part encode some "local elements of reality" . 

The non-degenerate case (n = 1) is generic and therefore the most frequently 
encountered in physical literature. Notably, U(l) Wilson loops are nothing else but 
the celebrated Berry's phase factor [2] accumulated by the wavefunction along a cyclic 
evolution. Since, in this singular Abelian case, covariance reduces to invariance. Berry's 
phase is potentially measurable. As a matter of fact, it has been directly observed 
in a variety of interference experiments [10-14]. This phase has demonstrated its 
usefulness to interpret some concepts or phenomena [15] such as, for instance, the 
paradigmatic Aharonov-Bohm effect [16] whose description in terms of Berry's phase 
was made by Berry himself. We can also mention the macroscopic polarization in 
crystalline dielectrics [17-19] expressed through Berry's phases of Bloch states across 
entire Brillouin zones [20] (also called Zak's phases), or the pumped charge in Cooper 
pairs pumps [21-23]. Moreover, locally. Berry's curvature field may be taken into 



I Wilson loops are often defined as the trace of holonomies. We choose an alternative definition, used 
in several contexts: Wilson loops as matrix representations of holonomies. 
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account to construe some physical effects. It is the case of the anomalous velocity 
of Bloch electrons, a correction to the usual term — coming from the band energy 
dispersion 8{k) — which manifests itself when the crystal momentum k is moved by 
external forces [24-27]. 

There is a close relation between Berry's connection and current operators given by 
partial derivatives (or gradient) of a Hamiltonian H with respect to parameter (s) [28]. 
This one is at the root of all the great geometrical and topological properties of systems 
endowed by such a current observable. Physical realizations are potentially multiple. 
For example, the current of Bloch electrons is proportional to dk,H , while the current 
crossing a phase polarized Cooper pair pump is proportional to d^pH , where is the 
superconducting phase bias of the device. It can be shown that, when we consider 
a non-degenerate eigenstate of a current operator of this form splits into a usual 
dynamical term (coming from the energy) and a geometric correction which can be in- 
terpreted in terms of Berry's curvature (see the example of the anomalous velocity). 
To the latter contribution is assigned the generic name adiahatic pumped current. In 
the context of the famous Thouless pumping [29], averaging the pumped current over 
the whole parameter manifold (a 2-torus) leads to a quantization of the charge. This 
stems from the topology of the line bundle over the torus, characterized by an invariant 
integer called first Chern number [1,6,28]. It is, mutatis mutandis, the same topological 
quantization as the Hall conductance for the integer quantum Hall effect [30,31]. Obvi- 
ously, that kind of topological invariance is highly interesting for metrological issues [32]. 

A. Joye et al [33] studied recently the current pumped in the degenerate non- 
Abelian case, too. They exhibited the differences between the Abelian and non-Abelian 
cases in a rigorous and somewhat abstract analytical approach. In this paper, we 
adopt the geometric counterpart of their work. Our aim is to give a more simple and 
comprehensive picture of adiabatic pumping, mainly for practical computations and 
later hypothetic applications (some predictions were already done within the framework 
of superconducting circuits [34,35]). We reserve section 2 to a review of the basic 
gauge principles underlying adiabatic dynamics, and we will emphasize the central role 
played by Wilson loops, especially when discretized for numerical purposes. Then, 
section 3 is devoted to the main subject. Considering a general n-fold degenerate level, 
we will express any pumped charge in terms of some Wilson loops' element. Contrary 
to what Wilson loop may seem, cyclic evolutions of the parameters are not required. 
The geometric difference between the two cases n = 1 and n > 1 will be identified. In 
the non-Abelian case, our holonomic formulation will allow a simple measure of "non- 
commutativity" of pumping cycles. Finally, in section 4, we make use of a "toy model" 
already encountered in non-Abelian pumping contexts [33, 34] to compare our previous 
results with simple analytical expressions. 
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2. Geometrical structure underlying quantum adiabaticity 

In this section, we briefly review the peculiar geometry of quantum systems adiabatically 
driven by classical parameters. Broadly speaking, such parameters can be degrees of 
freedom having their own dynamics (e.g. nuclei positions in the Born-Oppenheimer 
approximation), constraints eventually tuned by experimentalists (e.g. magnetic fields 
and voltages in the context of quantum circuits [36]), time t itself, etc. They appear 
as real arguments of the Hamiltonian, inducing a parameterized spectrum. If there is a 
total of m parameters x^, it is natural to represent the "parameter state" by a vector 
X = (x^, . . . , x™) in M™. Nevertheless, for example, a periodicity of H in can motivate 
us to compactify the "/i-th dimension" onto S^. Or, regarding some properties, defects 
may occur in a certain domain of M"^ inducing "holes" in it (e.g. gauge anomalies). As 
a generic consequence, we are likely to use an m-dimensional (real) smooth parameter 
manifold M, not necessary R™. Setting as the Hilbert space, we thereby construct 
a map H : M ^ ^{'H) assigning to p the Hamiltonian H{p). An evolution of the 
parameters (i.e. a path in M) induces an evolution of the wavefunction \1/ (i.e. a path 
in Ti) via the Schrodinger equation. In the well-known adiabatic assumption, under 
some hypothesis, the latter obeys to a simple parallel transport rule in an "eigenspace 
bundle" , as reviewed below. 

2.1. General settings 

Let us assume that the gap hypothesis is fulfilled for an isolated eigenvalue S of H, 
smoothly defined everywhere on M. That is to say, £^ is a smooth function on M and 
its finite multiplicity n is an invariant. Let (Bp ~ C" be the eigenspace corresponding 
to S{p) for each point p in M. In the adiabatic limit, if the state is initially found 
in CB^(o), it is well-known [37,38] that along a path 7 : r i— ?■ 7(t), the wavefunction 
"lives" at each time r in £-y(T-). The whole Hilbert space being useless, one reduces 
the structure when "gluing" together all the spaces (Bp to yield a smooth vector bundle 
TT : € — M with fibres vr^^(j>) = (Bp. In this bundle point of view, the adiabatic 
wavefunction \[^(r) G 7r~^[7(r)] moves horizontally to 7. More precisely, up to an 
overall phase factor, \I' obeys to the BWZ parallel transport rule built as follows. Let 
{\ua{p)) G 'n'~^{p)}a=i,...,n be a local smooth choice of orthonormal frames. Within 
this gauge fixing, the BWZ connection u is locally represented by a u(?T,)-valued 1- 
form field A: using a coordinate chart x = {x^, . . . ,x"^) = x^e^, A reads A^^dx^ with 



We now consider a path 7 : r 7(r) in M. In the adiabatic limit, the wavefunction 
is \I/(t) = e*''*^'^^ "^(t), where 



is the dynamical phase, and ip the parallel transport of ^1/(0) = ip{0) along 7. In 
components, if {ip) '■= ip°'\ua), the parallel transport equation is 
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Figure 1. Illustration of the gauge structure over a parameter manifold Af , in the 
Abelian case n = 1. Parallel transport in fibres along a loop C = 7t.o produces a 
holonomy sending ^/'(O) to ip{T). 



Its formal solution is thus 

where 7^^o = 7l[o,r] is the restriction of 7 to [0, r], and W is the Wilson operator acting 
on any path a through 

W[a]:=Ve~^^^ e\]{n). (1) 

W[a] is commonly called a Wilson line, or, if a is a (based) loop, a Wilson loop [39]. 
Finally, one obtains the standard expression of the adiabatic wavefunction [7]: 

|vl;(r)) = e^^^^^W\Ylrfl] /(O) KW)]). (2) 

Let us recall how A and W behave under a gauge transformation (i.e. a change of 
reference frames) generated by the smooth map g '■ p ^ g{p) G U(?t,): 

\Mp)) ^ \ua{p)) = \ub{p)) 9\{p)- (3) 
Each component transforms under (3) in accordance with the compatibility condition: 



Mp) ^ Mp) = 9 \p)Mp)9{p) + 9 \p) d^g{p), 
i.e. A (-T- g^^Ag + g^^dg. As for the Wilson operator, using the equality 

A^[7(r)]x''(r) = 1^[7,,o] 5.1^-^7^,0], 
it is easy to check that iy[7^Q] transforms according to the similarity rule: 



(4) 



(5) 
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Notably, if 7^,0 is a loop C, i.e. if 7(T) = 7(0) =: po; W[C] transforms covariantly with 
respect to its base point po '■ 



Properly, in the basis {\ua{j)o))}, W[C] represents the parallel transport map along C, 
itself an element of the BWZ holonomy group based at p^. In particular, when n = 1, 
W[C] reduces to Berry's phase factor which actually depends only on the unbased loop 
corresponding to C (see figure 1). Although it is abusive, by analogy with the Abelian 
case, W\C] is often called a non- Abelian geometric phase when n > 1. Of course, 
physical relevance is assigned to its eigenvalues which are gauge invariants of the parallel 
transport map. A good strategy consists in expressing measurable quantities as their 
functions (see [40] in the framework of quantum pumping). In the next section, however, 
it will be possible to link the pumped current to a single diagonal element of some Wilson 
loops. 

We have not yet spoken about the BWZ curvature Vt. Within a gauge fixing, it is 
locally represented by the 2- form field F := dA + ^ A ^ = | F^^ dx^ A dx" , with 



i.e. T ^ g Tg. 

2.2. (Discrete) Wilson loops as fundamental ingredients of the theory 

Because curvature fields and holonomies behave like observables, both are susceptible 
to carry some "elements of reality" , as we shall see in a while, within the framework 
of the pumping process. However, we will regard the latter as the most fundamental 
objects for two reasons. 

(i) Any component of the curvature tensor can be expressed by means of an 
infinitesimal Wilson loop. Indeed, in coordinates, consider a small rectangular 
loop c based at x and oriented along the /x-th then the i/-th directions. Call £^ its 
length in the /i-th direction and ^y in the iz-ih one, with £^ ~ (see figure 2). In 
the appendix, we derive the following relation: 



where the O symbol must be here understood in a matrix norm sense, and ^ ~ 
stands for the order of the "plaquette" size. The above formula can be reduced from 
the anti-Hermitian nature of F,.y: 



W[C]^W[C]=g^\po)W[C]g{po). 



Ffj^u — d^Ay — dpA^ + \A^, Ay\ 
the antisymmetric curvature tensor. It is obviously gauge covariant: 
F^,u ^ F^y = g-^F^yg = kdg-i{F^,^), 




(6) 




(7) 
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Xz = X + IvEv X2 

t * T 

1 >- • 

xo=x xi=x + ^^e^ 

Figure 2. The value of the curvature componcut Ff^^{x) is contained in Wilson loops 
of small oriented rectangles c, as defined in text. The proof given in the appendix 
needs a decomposition into four oriented segments £kj- 



(ii) A numerical scheme can easily be implemented to compute Wilson loops arising in 
our context. First, note that the naive discretization of formula (1): 

N 

W[C] = hm pTT (l - A^{p,)Ax';), 

with po = pn and Axj := Xj—Xj_^, seems inapplicable since it underlies smoothness 
in the gauge fixing, while numerical diagonalization routines "choose" eigenvectors 
in a somewhat erratic manner. This difficulty is overcome by the contextual 
definition of A^, leading to 



J b 

After introducing the overlap matrices S{p,q), whose elements are S"'^{p,q) := 
{ua{p)\uh{q)) , one obtains a formulation of great value for numerical treatments: 
the discrete Wilson loop 

N 

W[C]=\im vX[S{p,,p,.^). (8) 



In the special case n = 1, the above formula reduces to the discrete counterpart 
of Berry's phase factor: it can be seen as asymptotic versions of the Pancharatnam 
phase factor [41] or the Bargmann invariant [42]. Covariance of (8) is easily checked. 
Furthermore, it allows erratic frame choices along the loop: they cancel out in pairs, 
whereas the geometric information is solely extracted. In fact, one has removed the 
rigid constraints of differentiability and continuity in frame choices, properties which 
in fine reveal themselves unnecessary. Building an efficient algorithm based on formula 
(8) is a simple task: it needs basically scalar and matrix product evaluations, and 
a diagonalization routine for the Hamiltonian [43]. For = 1, it has been already 
implemented in miscellaneous contexts [17,18,23,44,45]. For n > 1, the result depends 
on the basis {\ua[po))} chosen on input. Then, in both cases, we let the program 
"choose" the frame {\ua{p))} of H{p) along the loop. We mainly apply formula (8) to 
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compute W[c\ in (7) as 

4 

W[c]=vX[S{x,,x,^r) + 0{l), (9) 
i=i 

where xq = Xi = x, xi, X2 and xs are c's vertices (see figure 2). 
3. Adiabatic charge pumping: a geometric viewpoint 

In this section, we exhibit the geometrical significance of adiabatic charge pumping. 
Such a process needs a quantum system equipped with (i) a current operator X; (ii) 
a Hamiltonian dependence in a controllable "pumping parameter" A G /, where / is 
an interval of M and (iii) a suitable Hamiltonian gauge such that X reads j^d\H [28], 
where g is a unit of charge. Physics of pumping becomes richer in the presence of 
other parameters. That is why we will assume H defined over an (m + 1) -dimensional 
parameter manifold M x I. With the current operator is naturally associated a mean 
transferred charge Q through the relation Q = (X), that is (in units of q) 

Q{r) ■.= \{nr)\dxH[p{r)-m{r)). 
3.1. Analytic derivation of the pumped charge 

For our purpose, we keep the pumping parameter A constant throughout quantum 
evolutions. Once given a path 7 in M, we define the set {7'^-'} of paths living in 
"horizontal" subregions A = est by: 7*^'(r) := (7(t); A). Along these, the wavefunction 
\l/ is well defined as a differentiable function of r and A. Consequently, we can write 

{^\dxH\^) = dx{^\H\^) - {dx^\H\^) - {^\H\dx^) 

= dx{H)-2hQ{{dr^\dx^)}, (10) 

where use has been made of the Schrodinger equation. We now assume an adiabatic 
evolution for a wavefunction belonging to an n-fold degenerate level. It is convenient 
to choose an initial frame such that ip°-[0]X) = 5"'', with 1 ^ r ^ n, and to set 
At-{t; \) := A^['j'-^''{t)]x^{t). Inserting the adiabatic approximation (2) in (10) and 
using the two relations below: 

drWH:^] = - A^{r;\)W[^l% 

5aW^[7|^]= - rW-'[^%]dxAr{T'-,X)W[^%]dT', 
Jo 

one obtains, after a little algebra, the instantaneous dynamical and geometrical (or 
pumped) contributions to the transferred charge. They are respectively 

dgd(r) =^dx£b''\r)]dT and (11a) 



dQr(r 



W~'H:^o]Frx{r;\)W[j^:^^,] dr, {lib) 



with F^x ■■= drAx - dxA^ + [A^, A 
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The dynamical term (11a) is obviously gauge invariant, local, and depends on the 
time parametrization (i.e. the dynamics of the path). However, it is independent of the 
initial state in the degenerate subspace. We now turn our interest to the geometrical 
contribution (11^). It is clearly sensitive to the initial state (r-label) and relies on an 
effective gauge structure in two variables r and A. We also recognize a component of 
the so-called (effective) twisted curvature tensor ^ [46]. Namely, 

5.,(r; A) := W~'[^^:^^,] F^,{t; A) W[^^:^^,]. (12) 

Under a gauge transformation g : (r; A) i— g{r; A), this one behaves as a Wilson loop 
based at (0; A): 

drxir; A) ^ §,A(r; A) = g-\0; A) drx{r; A) ^^(0; A). 

But, since one fixed ip°'{0] A) = 6°'^', the initial gauge freedom is constrained by 
g^'a{0; A) = g°'j.{0] A) = 6a,r for each index a. That makes the right-hand side of 
(116) gauge invariant. Second, anti-Hermiticity of Fr\ (and so of ^r\) guarantees 
the reality of Qg^- Third, its geometrical nature is due to its invariance under any 
time reparameterization r H- o"(r). Indeed, Wilson loops are themselves insensitive to 
parametrization, and one has 

dx^ d(T dx^ 
F^x dcr = Ff^x -r- 3- dr = F^a dr = Frx dr. 
do" dr dr 

In the Abelian case, ^ = F and dQg' is local, otherwise it is path dependent: the 
instantaneous pumped charge at any time is related to the whole history of the 
wavefunction. We point out that results (11a) and {lib) are easily generalizable for 
a vectorial current I = d\H with A = (Ai, A2, . . .). 



3.2. Geometric interpretation of the pumped charge 

In the previous paragraph we showed that the instantaneous pumped charge between r 
and r + dr is purely geometric and relies on ^, a tensor which behaves under a gauge 
transformation as a Wilson loop based at (0;A). These observations motivate us to 
search for a holonomic interpretation of Qg^- 

Let us first have a look at the charge pumped during an arbitrary small time interval 

At: 

AQ[^\T) = trrXrir;X)Ar + 0[iArn 

We may take advantage of formula (6) to compute the component -FT-A(r; A) in (12). 
If we denote by c the small rectangular loop pictured in figure 3, based at (r; A) with 
At ~ AA ~ A, equality (12) becomes: 

drxir; A) = W^-MtI^]^^^ W^bl^] + 0{A') 

= ^{a-w^ra-'°co7|;;)]}+0(A3). 
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Figure 3. (a) The parameter manifold M is symbolized by a horizontal plane while 
A's axis is seen in the vertical direction. The relevant paths 7'^' are made at constant 
A. It is shown in the text that the instantaneous pumped charge can be extracted from 
a loop (7x'^o)~^o c o "f^^Q, where c is a small rectangular loop oriented along the path 
and A's axis, (b) The effective space in variables r and A. 



Setting C the loop (7^^o) ^oco^^^q based at (0; A), one obtains an alternative formulation 
of the pumped charge: 

Ag«(r) = ^s{w^;[C]}+0(A2). (13) 

Note that in the Abelian case W\C] is equal to W[c\ and the rectangle c alone suffices 
for the computation. That represents geometrically the difference between the Abelian 
and no n- Abelian pumpings, in terms of locality or non locality. Formula (13), when 
combined with (8) and (9), is very useful for numerical purposes, as we will see later 
with an explicit example. 

As a second step, we would like to derive a similar expression for finite Ar. To this 
end, we will consider the charge Qg[7|^o^p] pumped along 7'^°' between any two times tq 
and Ti: 

QVbrS^=' r rir; \o) dr. (14) 

J TO 

Let {(x io I (''"' ^ b~o^ '^i] ^ ^} be the set of "vertical" segments oriented from (r; Aq) to 
(r; A) in the effective space (see figure 4). Then, we define rectangular loops c(r; A) as 
follows: 

c(r;A):=(Clt)"^o(7j^J-ioCl:l„o7(^). 
Setting T(r; A) := Vr[7|^„ o ([y^], the Wilson loop W[c{t; A)] can be expressed as 

W[c{t; A)] = T-\t; A) W[Ci]l] ^(^^ ^o), (15) 
where we have noticed the property 

T(r;Ao) = W^[7^^„']. (16) 
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Figure 4. An illustration of the various oriented segments and loops defined in the 
text to compute the pumped charge between times tq and ti. 



Thanks to the similarity rule (5), Vr[c(r; A)] in equation (15) appears to be the image of 
^[Ca^Ao] under the gauge transformation generated by T [47]. As for it transforms 
according to the compatibility condition (4): 

It is straightforward to verify that the partial derivative of A\ with respect to r is merely 
equal to F^-x = T~^FrxT. Or, in the integral form: 

Mr; A) = r T-i(r'; A)F,,(r'; A)T(r'; A) dr'. (17) 

J To 

The last result allows a curvature formulation of the Wilson loop iy[c(r; A)], through a 
"A-ordered" surface integral: 

W[cir; A)] = W[C,io] = e"^'o ^^(-^^')<^^' 



Vxexp{-1 [ T-'F,xTdrAdx], 

I Jy;(t;X) ) 



:(r;A) 

where S(r; A) is the well-oriented rectangular surface suspending the loop c(r; A) in the 
effective space. This is en passant an apparition of the so-called non-Abelian Stokes 
theorem [46,47]. The connection between the pumped charge (14) and our loops is 
established when taking c(ri; A) for a small difference AA := A — Aq: 

W[c{n; A)] = p.e-^'r"^-^-^^^')'^^' = t - Ax{n; Ao) AA + 0[(AA)2]. 
Using equation (17) and property (16), one obtains 

W[c{n; A)] = l-AA/"V-i[7;^j]F,A(r; Ao)W^[7^^,!]dr + 0[(AA)2]. (18) 



We set C = (Tt-o'o) ° '^i'^u A) o Tt^^, whose corresponding Wilson loop, deduced from 
(18) and (12), is 

W[C] = U - AA Hdrxir; Aq) dr + 0[{AXf]. (19) 

J TO 

After comparing (19) with (14) and using again anti-Hermiticity of ^^t-a, one obtains the 
final expression: 



Q'^'Ki.] = ^^{w:.[c]} + oiAX). 



(20) 
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Figure 5. Relative couplings between the excited state |e) and the three ground states 
Iffi); I52) and 1^3), in the four-state model discussed in section 4. 



We achieved our purpose: computations of pumped charges need in any case a 
simple Wilson loop evaluation by way of relation (8). Let us conclude this section with 
a consideration on loops in parameter manifold: non-Abelian charge pumping will give 
rise to non-commuting loops in the sense that loops T^i 



and Fg'^' based at the same 

point will generically be such that Q'^;^[V^^'' o T^"] ^ Q(.'-'[r^'' o r^"']. If Q^fr'^'l and 



J g i - - -Z i / g 

Qg'[r2'*^] rely respectively on loops Ci and C2 in equation (20), one can measure their 
"non-commutativity" through 

1 



AA 



W[C,],W[C2 



0(AA). 



Some experiments are already proposed to demonstrate the non-commutativity of 
pumping cycles [34,40]. 



4. An illustrative example 

We illustrate the results derived in the previous section with a well-known model 
appearing in atomic [9] and superconducting [34] systems. It consists of an excited 
state |e) parametrically coupled to three independent ground states \gi), \g2) and Ig^,). 
In [34] the Hamiltonian reads 

H{r;\) = Ee\e){e\ +x|e)(^i| +y\e){g2\ +ze'^ \e){gs\ +h.c., 

where x, y and z are three real parameters, r is the vector (x, y, z) in the parameter 
space M"^, and is the energy assigned to state |e) (eventually zero), while energies 
of the ground states are set to zero (see figure 5). In the context of [34], H{r] A) is a 
truncated Hamiltonian in the Coulomb blockade regime where |e), \gi)i \g2) and \g^) are 
the four relevant charge states (i.e. Cooper pairs in excess on different coupled islands). 
The parameters are then Josephson couplings between islands, the pumping parameter 
A is the phase bias between two superconducting leads and the current operator is 
the Cooper pair current from one lead to the other. Experimentally speaking, the 
parameters x, y and z live only in a restricted range of values. One may ask r to belong 
to a connected subset of but this restriction is somewhat irrelevant. When r 7^ 0, 



The Geometry of (non)-Ahelian adiahatic pumping 13 



z 




Figure 6. The charge pumped along a path F in R'^ is equal to the one pumped along 
r's radial projection onto the unit sphere S^. Here, F is a loop and its projection P(r) 
a circle. 



the Hamiltonian has three distinct eigenvalues: Eq = and E^, the second being 
doubly-degenerate over M^\{0} while the non-degenerate eigenvalues E± depend only 



on r := ||r|| = ^^/x^~+^y^~+^ : 



2 

In the special case r = 0, the level Eq becomes at least three-fold degenerate: the 
origin is thus a defect for the U(2) gauge structure and we reduce the parameter space 
to R'^\{0}. In order to perform non-Abelian pumping we focus on the eigenspaces 
associated with Eq. A smooth choice of orthonormal frame |^2)} can be made 

over R^\(Oz) x / as follows: 

I / \x\ y \9l) - X \g2) fr,-, ^ 

\Mr; A)) = (21a) 

_|_ yZ 



-ze 



x\9i)+y\92)) + ix^ + y^)\93) 



\u2{r; A)) = ^7^^^ ^'^^^^ ' . {21b) 

y (x^ + y'^){x^ + 2/^ + z^) 

A similar and compatible choice can also be found over e.g. M^\(0?/) x / to cover entirely 
the parameter space: 

ze'^\gi) -x\gz) 



\ui{r;\)) 



I y^^-^ ^ -yjx \gi) + ze '^Igs)) + {x^ + z^)\g2) 

We note that all these states are invariant under a rescaling r H- kr, k > 0. 
Therefore, it is sufficient to retract ]R'^\{0} onto its embedded unit sphere S^: the 
pumped charge along a path 7 in will be equal to the pumped charge along its 
radial projection ^(7) onto (see figure 6). Using local spherical coordinates on S^: 
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r = (sin6' cos <y9, sin6' sin cos6'), the eigenstates (21a) and (216) read: 
\ui{r; A)) = sin Ifi-i) - cos \g2) 

\u2{r] A)) = — e^^ cos 9{ cos ip \gi) + simp \g2)) + sin 9 Ig^). 
Within this choice, the components of the effective gauge fields have simple expressions: 

-e'^ 







and: 



A,{t;X) = -cos[e{T)]ip{T) 
A,{t;\)=icos''[9{t)] 

sm[2e{r)]f sin[e(r)]0(r)e*^ 

j^tx{t; a) = 




2i \ sm[e{T)]ip{T)e~'^ 2^(r) 

The simplicity of the above relations allows an exact computation of the Wilson line 
involved by the path followed between times and r: 



where we set: 



«r,o = / cos[6'(r')] V9(r') dr' = / cos6'dv9. 

Jo J■y-rA^ 



cos(a;T-_o) — e*'^ sin(aT- q) 
"''^ sin(aT-^o) cos(q;t-,o) 



One has thereby the basic ingredients to calculate via equation [lib) the transferred 
charges Q^^' and Qg' along 7*^'. Namely, 

g^^'[7|^] = sm\ar,o)sm'[9{T)], (22a) 
QfH?o] = cos2(«.,o) sin2[^(r)] - sin^[e(0)]. (226) 
Note that they do not depend on A's magnitude; here, A is a pure pumping parameter 
with no observable effects. It is straightforward to check that, once the vertical Wilson 
line 



^(C1:^aa,a) = ^Aexp I - A^{r; A') dA' j 



1 

Q g-icos2[6»(r)]AA 



is obtained, formula (20) yields to the same results (22a) and {22b). In particular, if 
at time T one has r{T) = r{0), the path 7^^^ is a loop F and the two pumped charges 
above are opposite: g^fp] = siii^iaTfi) sin2[^(0)] = -Qf [F]. 

We conclude that section by a numerical test in the present non-Abelian example 
(for an Abelian case, see e.g. [23]). The algorithm is built on the discussion of subsection 
2.2 applied to formula (13). To check its validity, consider a mere situation: a circle 
whose axis is {Oz) and such that cos 6* = 1/4. Whatever A be, if the circle is followed 
at constant velocity ip from r = to r = 1, expressions (22a) and [22b) become: 

Q'^n^o] = sin' (^-) = -Qrra. 

In figure 7, the above expression together with the numerical result is plotted, and we 
find that they perfectly coincide. 
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r 

Figure 7. Analytical and numerical plots of the pumped charge in the circular example 
given in the text, for radius 1, cos6' — 1/4 and A = 0. The algorithm is based on 
formulas (8) and (13). 

5. Conclusion 

We have shown the relevance of Wilson loops to describe adiabatic charge pumping, 
mainly for practical purposes. In that pure geometric approach, one may particularly 
take advantage of formulas (8) and (13) to compute pumped charges in any problem. In 
the future, it should be interesting to investigate the geometry of this process in open 
systems, and to have a look at topological issues in the non-Abelian case. 
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Appendix. Curvature extraction from rectangles 

In this appendix, we give a proof of formula (6). First of all, we denote by xq = x, 
Xi, X2 and xs the vertices of the oriented rectangle c, as depicted in figure 2, such that 
xi = X + i^e^ and x^ = x + lyCy. Then, we decompose the loop into four oriented 
segments joining its vertices: £io, ^21, ^32 and £o3) with l^j going from Xj to x^ (see 
figure 2). It suffices to treat £10 and to apply the result to the others. This path can be 
parameterized as follows: 

£10 : [0, 1] 9 s ^ £io(s) = Xo + s{xi -Xq)=Xq + s E M™. 

Within a gauge choice, its corresponding Wilson fine reads: 
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Expanding the path-ordered exponential as a usual Dyson series, one picks the 
expression of W^[^io] up to the second order in (in a matrix norm sense): 

W[iio] = l + y^i-£,)' dsJ ds2---/ 'dsJA,[£^oisi)]A^[hois2)]---A^[£,ois,)]} 
fc^l Jo Jo Jo ■' 

= l-i, f\sA,[£,o{s)]+£l f'ds rds'A^[ho{s)]A^[ho{s')] + 0{il) 
Jo Jo Jo 

= 1 - A,{xo) + \ {aI{x,) - d,A,{x,)]£l + 0(£y or 

= 1 - A,{x,) + \ {aI{x,) + d,A,{x^)] ll + 0{tl). 
Therefore, the Wilson loop can be expressed as the product 
W[c] = W[£^^] ■ W[£^2] ■ W[£2i] ■ W[ho\ 

= {l + A,{xo) L + \ {aI{x^) + d,A,{xo) } li) 

■ [t + A^{x,) £^ + i [aI{x,) + d,A,{x,)]e^ 

■ (l - A,{xo) £, + \ {a%xo) - d,A,{x,)]£^^ + 0{£'), 

where £ stands for the order of the rectangle size: £ ~ £^,£v Finally, using Ay{xi) = 
A^.{xo) + £ud^A^{xo) + 0(£2) and v4^(x3) = A^{xo) + £^d^A^{xo) + 0{£^), one obtains 

W[c] = 1 - [d^A,{xo) - d,A,{xo) + [A,{xo),A,{xo)] } £/, + Oi^) 

= 1-F^,{x) £/, + 0{£^). 
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